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ABSTRACT 



We present new observations leading to an improved black hole mass esti- 
mate for the Seyfert 1 galaxy NGC 4593 as part of a reverberation-mapping 
campaign conducted at the MDM Observatory. Cross-correlation analysis of 
the H/3 emission-line light curve with the optical continuum light curve re- 
veals an emission-line time delay of Tcent = 3.73 ± 0.75 days. By combining 
this time delay with the H/5 line width, we derive a central black hole mass of 



^Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210; 
denney, bentz, dietrich, ghosh, peterson, pogge@astronomy.ohio-state.edu 

^School of Physics and Astronomy, University of St. Andrews, Fife, KY16 9SS, Scotland, UK; emcl4, 
kdhl@st-and.ac.uk 

■^Department of Astronomy, University of Michigan, Ann Arbor, MI 48109-1090; fogel, dor, mg- 
walker @umich.edu 

''Current address: Physics and Astronomy Department, 3270 Biomedical physical Sciences Building, 
Michigan State University, East Lansing, MI 48824; kuehncha@msu.edu 

^Institute of Astronomy, School of Science, University of Tokyo, 2-21-1 Osawa, Mitaka, Tokyo 181-0015, 
Japan; minezaki, yoshii@ioa.s.u-tokyo. ac.jp 

^Present address: National Research Council Canada, Herzberg Institute of Astrophysics, 5071 West 
Saanich Road, Victoria, BC V9E 2E7; Christopher. onken@nrc-cnrc.gc.ca 

^Crimean Astrophysical Observatory, p/o Nauchny, 98409 Crimea, Ukraine; sergeev, 
vpronik@crao . Crimea, ua 

^Isaak Newton Institute of Chile, Crimean Branch, Ukraine 

^Steward Observatory, The University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721; mvester- 
gaard@ as . arizona. edu 

^^Research Center for the Early Universe, School of Science, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, 
Tokyo 113-0033, Japan 



- 2 - 



Mbh — (9.8 ± 2.1) X IO^Mq, an improvement in precision of a factor of several 
over past results. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 

1. INTRODUCTION 

Over the past two decades, reverberation mapping (Blandford & McKee 1982; Peterson 
1993) has emerged as the most widely applicable method for directly measuring the masses 
of supermassive black holes (SMBH) in active galactic nuclei (AGNs). The appeal of this 
method arises from its physical simplicity, despite the complicated physics within the central 
regions of these galaxies. Moreover, unlike other methods for measuring the masses of 
SMBHs, it does not require high angular resolution. On the other hand, this method is 
challenging because it requires a large series of observations, well spaced in time and taken 
over a long duration (at least three times the length of the longest physical timescale of 
interest). 

The reverberation-mapping method depends on continuum flux variations from the cen- 
tral source, presumably an accretion disk around a SMBH. Continuum photons are absorbed 
and then re-emitted as line photons in the broad line-emitting region (BLR). Variations in 
the broad-line flux have a time lag relative to the continuum, corresponding to the light- 
travel time across the BLR. Through measurements of this time lag, the size of the region 
can be determined, which then leads, through virial arguments and a measurement of the 
broad-line width, to the black hole mass, Mbh- 

To date, black hole masses have been determined for nearly 40 AGNs using reverberation 
mapping (see the compilation by Peterson ct al. 2004). While this demonstrates how far the 
field has progressed, Peterson ct al. (2004) make it clear that substantial work is still needed 
to decrease uncertainties in time lag, and thus Mbh, measurements in order to better test and 
calibrate the various black hole mass-dependent scaling relationships, such as the Mbh — cr* 
relationship, where cr* is the bulge velocity dispersion. This need for better reverberation- 
based mass estimates is the driving force behind the campaign we describe below. 

One requirement for a successful reverberation mapping program is good temporal sam- 
pling (Horne et al. 2004). The major goal of the current observing program is to improve 
the sampling rate for reverberation-mapped AGNs that were undersampled in previous cam- 
paigns. The mass of the black hole in NGC 4593 was previously measured by reverberation 
mapping (Dietrich et al. 1994; Onken et al. 2003), but the H/? lag determination of Tcent = 3.1 
days was unfortunately much smaller than the average interval between observations of 15.8 
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days. This resulted in uncertainties (+7.1, —5.1 days) much larger than the time delay itself, 
and as such, were consistent with no lag at all. By revisiting this object with our program, 
we were able to achieve the sampling necessary to drive down the previous high uncertainties 
in the time lag and black hole mass. Here, wc present new lag and Mbh determinations from 
reverberation mapping of NGC 4593. The observational uncertainties from this campaign 
represent a factor of several improvement over previous results. 

2. OBSERVATIONS AND DATA ANALYSIS 

Observations of NGC 4593 were obtained as part of a large reverberation mapping cam- 
paign undertaken in early 2005 on the 1.3m McGraw-Hill telescope of the MDM Observatory 
on Kitt Peak in Arizona. Details beyond the descriptions below, for our spectral and pho- 
tometric observations and data reduction, are given by Bentz et al. (2006b). NGC 4593 
is one of three AGNs from which we were able to observe sufficient variability within the 
time frame of our program to warrant a reverberation analysis. It was first recognized as a 
host to a type 1 AGN by Lewis et al. (1978) as part of the Michigan-Tololo Curtis Schmidt 
survey for extragalactic emission-line objects (see MacAlpine et al. 1979). NGC 4593 is a 
type (R)SB(rs)b spiral galaxy at a redshift oi z — 0.0090. 

2.1. Spectroscopy 

MDM Observations: We obtained spectra of the central region of NGC 4593 at 
24 epochs between 2005 February and April with the Boiler and Chivens CCD spectro- 
graph on the 1.3-meter telescope at MDM Observatory, at a spatial scale of 0'.'75 per pixel. 

For our campaign, we used a grating of 350 grooves/mm. corresponding to a dispersion of 
1.33 A/pixel. The spectral coverage was from ~ 4300 - 5700 A, centered on H/5A4861 and 
the [O III] AA4959, 5007 lines, resulting in a spectral resolution of 7.6 A across this range. 
The slit width was set to 5'.'0 projected on the sky, with a position angle of 90°. Figure 
1 shows the mean and rms spectrum of NGC 4593, made from the complete set of MDM 
observations. 

CrAO Observations: We also acquired spectra of the nuclear region of NGC 4593 
from the Nasmith spectrograph with the Astro-550 580 x 520 pixel CCD (Berezin et al. 
1991) on the 2.6-m Shajn telescope of the Crimean Astrophysical Observatory (CrAO). For 
these observations a 3'/0 sht was utilized, with a 90° position angle. Spectral wavelength 
coverage for this data set was from ~ 4300 — 5600 A, with a dispersion of 2.0 A/pix and a 
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spectral resolution of 8.2 A. 

2.2. Photometry 

MAGNUM Observations: In addition to spectral observations, we also obtained V- 
band photometry from the 2.0-m Multicolor Active Galactic NUclei Monitoring (MAGNUM) 
telescope at the Haleakala Observatories in Hawaii, imaged with the multicolor imaging 
photometer (MIP) as described by Kobayashi et al. (1998a,b), Yoshii (2002), and Yoshii et al. 
(2003). Photometric reduction of NGC 4593 was similar to that described for other sources 
by Minezaki et al. (2004) and Suganuma et al. (2006), except the host-galaxy contribution to 
the flux within the aperture was not subtracted and the flltcr color term was not corrected 
because these photometric data were later scaled to the MDM and CrAO continuum light 
curve (as described below). Also, minor corrections (of order 0.01 mag or less) due to the 
seeing dependence of the host-galaxy flux were ignored. 

2.3. Light Curves 

Following reduction of the data, light curves were created for the subsequent cross- 
correlation analysis. The spectral fluxes of the MDM observations were calibrated based 
on the [O ill] A5007 line flux between observed frame wavelengths of 5022 — 5062 A in the 
mean spectrum. Individual spectra were scaled to this reference spectrum using software 
that employs a goodness of fit estimator (van Groningen & Wanders 1992). Because 
few MDM observations were taken under photometric conditions, our absolute [O ill] A5007 
flux calibrations are not reliable. Therefore, the flnal light curves are calibrated to the 
[Olll] A5007 flux given by Dietrich et al. (1994). Light curves from the CrAO spectra were 
made following a similar method, only the flnal light curve flux measurements were not 
calibrated to the Dietrich et al. (1994) value because this data set was later scaled to the 
MDM data set, as described below. Continuum and line flux measurements were obtained 
by first fitting the continuum on either side of the H/3 fine, between wavelength ranges 
4795 - 4815 A and 5120 - 5170 A in the observed frame. The H/3 fiux was calculated by 
integrating above the continuum over the wavelength range 4825 — 4963 A. The continuum 
fiux density was then taken to be the average over 5120 — 5170 A. Although at least some 
contamination from Fell emission is unavoidable when measuring the continuum in the 
optical, we chose this continuum region because it is the cleanest window, with respect to this 
contamination, within the spectral coverage of our observations. In addition, Vestergaard 
& Peterson (2005) demonstrate that although Fell emission is variable, the amplitude of 
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variability is generally low. We do not expect Fe ll contamination to significantly affect our 
measurements. 

In order to intercalibrate the data from various sources (i.e. place them on a consistent 
flux scale) , corrections to the measured line and continuum flux values are necessary because 
of systematic differences, mostly attributable to seeing and aperture effects (see Peterson 
et al. 1991). Given data sets from two sources, flux differences were determined by comparing 
all possible pairs of simultaneous observations between the two data sets. Here, we relaxed 
the meaning of 'simultaneous' to include pairs separated by at most 2.0 days so as to increase 
the number of pairs contributing to each comparison. However, this relaxed assumption does 
not produce significantly larger (i.e. less than 1.0%) uncertainties, implying that there is no 
evidence for intrinsic variability on such short timescales. Based on the average difference 
in flux between these closely spaced pairs, a single multiplicative point-source correction 
factor, (^crAO = 0.89, was applied to all points in both the emission-line and continuum light 
curves of the secondary data set (CrAO) to scale it to the primary set (MDM), following 
the methods of Peterson et al. (1991). Physically, this scaling factor accounts for differences 
in the amount of [O III] A5007 flux measured due to different aperture geometries and seeing 
between the two data sets. Seeing can also effect the internal calibrations of the continuum 
and B./3 emission-line flux measurements, based on F([0 ill] A5007), which is assumed to 
be constant. However, given the large slit widths used in the current MDM and CrAO 
observations (S'.'O and 3'/0, respectively) as well as the compact (no larger than 1'.'7) emission 
region of [Olll] in NGC 4593 (Schmitt et al. 2003), seeing is an insignificant effect in our 
internal calibrations. Further tests of the effect of seeing on our internal flux cahbrations 
will be addressed in §3. 

In addition to this multiplicative correction, an additive correction factor, GcrAO — 
—1.60 X 10~^^ erg s~^ cm~^ A~^, was also made to the continuum flux, accounting for the 
starlight contribution in the different apertures of the various data sets. These correction 
methods were employed to scale both the continuum and H/9 light curves from the CrAO 
observations to the MDM data, since the MDM data formed the largest, individual set, and 
thus served as the base onto which to tie the others. For the same reason, the photometric, 
V-band observations from the MAGNUM data were then scaled to the newly combined MDM 
and CrAO continuum fight curve, with additive correction factor Gmagnum — 0.65 x 10~^^ 
erg s~^ cm~^ A~^. The continuum light curve was finally corrected for host galaxy starlight, 
where, for this slit geometry, the contribution is measured to be Fgai(5100A)— (Lllloiii) x 
10~^^ erg s~^ cm~^ from a recent Hubble Space Telescope image taken with the Advanced 
Camera for Surveys High Resolution Camera. The methods by which this flux contribution 
was determined are described by Bentz et al. (2006a), and the complete results specific to 
this object will appear in a future paper. The full set of observations used in constructing 
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these light curves can be found in Table 1, and Figure 2 shows the corresponding hght curves. 

3. TIME SERIES ANALYSIS 

For the time series analysis, the light curves from Figure 2 were modified by first merging 
the three sets of observations by binning data points more closely spaced than 0.5 day 
through the use of a variance-weighted average. Cross-correlation analysis was restricted to 
the subset of the light curves where both continuum and emission-line measurements are 
available, i.e. JD2453430 - JD2453472. In addition, two observations from the CrAO data 
set, JD2453463.4 and JD2453465.9, show anomalously high continuum flux for reasons that 
could not be determined. Conducting the cross-correlation analysis including these points 
resulted in lag determinations consistent with removing the points, although the uncertainties 
in the lags were higher. Since this time period was otherwise well sampled, none of the other 
observations show a similar peak in continuum flux, and because the results from exclusion 
are consistent with inclusion, we omit these two continuum points from the final light curve. 

Figure 3 shows the modified light curves that were used for the cross correlation analysis. 
The statistical parameters describing these final light curves are shown in Table 2. Column 
(1) gives the spectral feature, and Column (2) shows the number of data points in each hght 
curve. Columns (3) and (4) are mean and median sampling intervals, respectively, between 
data points. The mean flux with standard deviation is given in column (5), while column (6) 
shows the mean fractional error, based on comparison between closely spaced observations. 
Column (7) gives the excess variance, calculated as 

where cr^ is the variance of the observed fluxes, 5^ is their mean square uncertainty, and 
(/) is the mean of the observed fluxes. Finally, column (8) is the ratio of the maximum to 
minimum flux in the light curves. 

We conducted the time-series analysis by cross correlating the continuum light curve 
with the H/3 light curve using two methods designed for unevenly spaced observations. The 
primary method uses an interpolation scheme that averages two results: flrst, cross cor- 
relating an interpolated continuum light curve with the original emission-line light curve, 
and second, cross correlating the original continuum with an interpolated emission-line light 
curve (Gaskell & Peterson 1987). We employed an interpolation interval of 0.5 day (equal to 
half of the median time span between observations) , chosen to create equal spacing on a scale 
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smaller than that of the actual data, yet large enough to prevent the introduction of artifacts 
due to the correlation of adjacent interpolated flux values. The correlation coefficient, r, is 
computed from pairs of points from each light curve matched based on different lags being 
imposed upon the trailing curve. For a given time lag realization, unpaired points at the 
beginning or end of the time series are excluded from the correlation analysis. By calcu- 
lating r for multiple potential lag times, a cross-correlation function (CCF) is constructed, 
shown in Figure 4, which gives the correlation coefficient at each of these different lag times, 
r. From the CCF, we can characterize the time delay through two parameters. The peak 
lag, Tpeak, is simply taken as the lag time which produces the highest correlation coefficient, 
''max, whereas the centroid lag, Tcent, is the centroid of the CCF, computed from values with 

r > O.Srmax- 

The uncertainties in the time lag measurements were computed via model-independent 
Monte-Carlo simulations and based on the bootstrap method referred to as FR/RSS (for 
flux redistribution/random subset selection), described by Peterson et al. (1998) and which 
utilized the modifications of Peterson et al. (2004). A large number of realizations of this 
method (e.g. 10,000 in this study) are employed to build up two distributions of CCFs: the 
cross correlation centroid distribution (CCCD) and the cross correlation peak distribution 
(CCPD). The values Tcent and Tpeak are the means of these respective distributions. The 
uncertainties are calculated such that 15.87% of the reahzations yield values larger than the 
mean plus the upper error, and 15.87% yield values smaller than the mean minus the lower 
error; these would correspond to the la errors in the case of a Gaussian distribution. Table 3 
gives the calculated values of Tcent and Tpeak for NGC 4593, after being corrected for time 
dilation. 

The second cross correlation method we employ creates a discrete correlation function 
(DCF), described by Edelson & Krohk (1988), including modiflcations by White & Peterson 
(1994). Shown in Figure 4, the DCF is created by determining correlations at different 
lag times between the H/3 and continuum light curves, just as in the interpolation method, 
but through time binning of data rather than interpolation. Under-sampled data or data 
with large time gaps could lead to spurious lag determinations through the interpolation 
method, in which case there are advantages to using the DCF method. For example, only 
the actual data are used in the analysis and due to the binning, statistical uncertainties 
can be assigned to the correlation coefficient for each bin (see White & Peterson 1994). 
For this analysis, a time bin of 2.0 days was adopted. One can see from Figure 4 that the 
DCF agrees well with the CCF from the interpolation method showing that interpolation is 
not introducing artifacts. Figure 4 also shows the auto-correlation function (ACF), which is 
computed by cross correlating the continuum with itself. The ACF is interesting because the 
CCF is the convolution of this distribution (ACF) and the transfer function, which describes 
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the emission-line response to the continuum variations and is the quantity reverberation 
mapping is fundamentally attempting to recover (see Peterson 1993). 

The cross-correlation functions in Figure 4 also demonstrate that our continuum and 
H/3 emission-hne internal flux calibration to the [O ill] A5007 emission hne are not adversely 
affected by nightly changes in seeing. Because calibration errors due to seeing would result in 
correlated errors, with both the contimiiim and the line flux affected in the same direction, 
if variable seeing is a significant effect, it would lead to a spike in the cross-correlation 
function at a lag of zero days. Figure 4 clearly shows no such feature in either the CCF 
or DCF for NGC 4593. Bentz et al. (2006b) report on results for NGC 4151 from this 
same observing campaign and do not see any evidence for this zero lag spike either. The 
only object from this campaign in which a zero- lag spike is detected in the cross-correlation 
function is NGC 5548 (Bentz et al. 2006c), which was in a very low-flux state at the time. 
In this case, the amplitude of variability was so low that the weak correlated-error signature 
was above the detection threshold. 

To investigate the possibihty of detecting bulk motions within the BLR, we then broke 
the H/? broad emission-line into velocity bins and performed the cross correlation anaylsis 
between these individual bins. First, the H/3 line was divided in half based on the peak in the 

mean spectrum (see Fig. 1), with the blue side defined by wavelengths 4825 - 4903 A, and 
the red side defined by wavelengths 4903 - 4963 A. The fiux was determined by integrating 
the line fiux above the same continuum fit as described above. The lags determined through 
this cross correlation analysis are Tcent = 0.5 ±1.0 day and Tpcak = O-^^i's day, with variations 
in the blue side slightly lagging behind those in the red side. Second, the line was broken into 
a 'core' and 'wings' such that the integrated flux in the core was equal to the summed flux 
in both red and blue wings. These regions were deflned from the mean spectrum such that 
the blue wing constituted the integrated flux between wavelengths 4825 - 4881.8 A, the core 
between 4881.8 - 4924.2, and the red wing between 4924.2 - 4963 A. The lags determined for 
this scenario are Tcent — 0.5 ± 0.7 day and Tpcak = 0.1 ± 1.0 day, with variations in the wings 
leading the core as would be expected. The CCFs from both of these analyses are shown 
in Figure 5. Unfortunately, any signature of infiow/outfiow or bulk Keplarian motions was 
too weak to be detected with confidence through this investigation, as the lags measured for 
both scenarios were consistent with zero. 



4. BLACK HOLE MASS 



Assuming that the motions within the line-emitting region are gravitationally dominated 
so that the virial theorem can be utilized, the mass of the black hole can be defined such 
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-'WBH — ^ 

where r is the measured emission-hne time delay and AV is the emission-hne width, which 
can be characterized by either the full width at half maximum (FWHM) of the broad 
emission-line or by the emission-line dispersion, crune- The dimensionless factor / depends on 
the structure, kinematics, and inclination of the BLR and is of order unity. By normalizing 
the reverberation-based black hole masses to the Mbh — o"* relationship of quiescent galaxies, 
Onken et al. (2004) have shown that / has an average value of 5.5, when the emission-line 
dispersion, or second moment of the profile, is used (as opposed to the FWHM^) for AV^. 

After first removing the narrow line component of the H/5 emission line, the FWHM and 
the line dispersion were measured from both the mean and the rms spectra of NGC 4593. 
Peterson et al. (2004) describe in detail the bootstrap method used to determine these 
quantities and their respective uncertainties. The measured values for the FWHM and crune 
are given in Table 3. 

To calculate the black hole mass for NGC 4593, we use the centroid lag, Tccnt, for 
the time delay, r, and the line dispersion, crune, of the B./3 emission line, measured from 
the rms spectrum, for the emission-hne width, AV. Peterson et al. (2004) argue that this 
combination of parameters gives the most reliable black hole mass determinations, based on 
virial arguments and fits between all possible combinations of parameters. The black hole 
mass calculated for NGC 4593 from this work is Mbh = (9.8 ± 2.1) x IO^Mq, where the 
uncertainties quoted include statistical and observational considerations but not an intrinsic 
uncerainty in the method, no larger than a factor of 2-3, that accounts for unknowns such 
as inclination. 

Our result is consistent with the mass estimate based on a previous reverberation pro- 
gram. Onken et al. (2003) reanalyzed observations of NGC 4593 by the "Lovers of Active 
Galaxies (LAG)" consortium (Dietrich et al. 1994) and determined a black hole mass based 
on H/3 to be Mbh = i^Atj]"^) x lO^M©, where we have used the Onken et al. (2004) calibra- 
tion (Peterson et al. 2004). The improved time sampling of the new reverberation mapping 
observations has significantly decreased the uncertainties in the black hole mass estimate for 
this object. 



^See Collin et al. (2006) for a discussion of systematic differences between the use of the FWHM and the 
line dispersion in calculating the virial product, and thus Mbh- 



(2) 
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5. DISCUSSION AND CONCLUSION 

Results presented here for the emission-hne lag and black hole mass of NGC 4593 of 
Tccnt = 3.73 ± 0.75 days and Mbh = (9.8 ± 2.1) x IO^Mq, for these respective quantities, 
represent a factor of several improvement over past measurements. Previous measurements 
for this object from LAG data exhibited much larger uncertainties due in part to the average 
sampling interval, which due mostly to bad weather, was much larger (by about a factor of 4) 
than the emission-line lag of this object. In addition, for rather short observing campaigns, 
some of the success of reverberation mapping comes from serendipitously observing a large 
variability event during the campaign. The best results come from observing a large increase 
and then decrease (or vice versa) in the flux, which we observed in this campaign but was 
unfortunately not seen to the same degree in the previous campaign of NGC 4593. This 
represents yet another example exhibiting the need for longer observing campaigns. 

The results reported here were obtained as part of a larger spectroscopic monitoring 
campaign whose primary goal was to improve the emission-line lag and black hole mass 
measurements for some of the nearest, apparently brightest AGNs. The proximity of these 
AGNs makes them especially important not only as calibrators of the AGN Mbh^ct* and BLR 
radius-luminosity relationships (Bentz et al. 2006a), but also as candidates for measurement 
of their black hole masses through other means that depend on high angular resolution. Key 
elements of this program were (1) scheduled daily observations of each of these relatively low- 
luminosity AGNs, (2) supporting observations at multiple sites to mitigate the effects of gaps 
in the time series caused by weather, and (3) high-quality homogeneous data that permitted 
relative spectrophotometric flux calibration at better than the 2% level. Our 42-night pro- 
gram yielded improved HP lags and black hole masses with a factor of several improvement 
in precision relative to previous campaigns for three of the six AGNs in our sample, NGC 
4151 (Bentz et al. 2006b), NGC 5548, which was observed in the lowest-luminosity state yet 
recorded (Bentz et al., in preparation), and NGC 4593, as described here. Three other AGNs 
in our program, NGC 3227, NGC 3516, and NGC 4051 were insufficiently variable during 
this relatively brief campaign, underscoring the point that longer duration programs are 
necessary for detection of variability signatures that are favorable for reverberation analysis. 
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Table 1. Continuum and H/3 Fluxes for NGC 4593 



JD 
(-2450000) 


Fa (5100 A) 
(10~^^ erg cm~^ A~^) 


H/? A4861 
(10~^^ erg cm~^) 


Data Set 


ooyi.yoD 


10.69 ± 0.10 




iViA*jiN U iVi 


o4ii.yzD 


9.23 ± 0.10 




i\/r A r^i\TTTi\/r 

iVlAljiN U iVi 




9.86 ± 0.05 




iViriVjiN U iVi 


■^/i "^n 8/1/1 


8.94 ± 0.18 


/I KQ _i_ n no 
4.oy zt u.uy 


iViUiVi 


Q/l QD 


9.78 


± 0.08 




iVl/\<ori\ U iVi 


040i .OOZ 


8.75 ± 0.18 


4.0Z It U.uy 




O400.oi0 


8.52 ± 0.17 


4.oi ic u.uy 


A/T'niv/r 

iViUiVl 


Q/l /I QQ 

o4o ( .4yy 


9.68 ± 0.32 


A R7 -1- n 1 n 

4.0 ( It U. iU 


Pt- a n 


o4o ( .yzo 


9.39 ± 0.19 


A 77 4- n 1 n 

4. ( ( It U. iU 


A/rF>A;T 
iViiJiVi 


Q/l Q8 07^ 


10.37 ± 0.04 




iVi/i*jiN U iVi 


■^/l '?8 80"^ 
O4OO.OU0 


9.81 


± 0.20 


O.UO It U.IU 


iViUiVi 


"^A "^Q 78Q 

o4oy. / oy 


9.60 


± 0.19 


O.UU It U.IU 


iViUlVi 


OA Art QOA 
o44U.oZ4 


9.51 


± 0.19 


O.Uo it U.iU 


iVlUiVl 


Q/1/11 8/18 
044i .o4o 


9.24 


± 0.19 


c; 1 9 _i_ n 1 n 

0. iZ It U.iU 


i\/rF»A/r 


Q/l/IO 8on 
044Z.oZU 


8.92 


± 0.18 


K OA -I. n 1 1 
O.o4 It U. i i 


iViUiVi 


Q/l/l'? 890 
0440.oZU 


9.05 


± 0.18 


K O'i -4-0 11 
O.Zo It U. i i 


iVlUiVl 


OA A A A en 
C)444.4oU 


8.75 


± 0.29 


K n8 1 nil 
O.Uo it U.il 




O440.40D 


7.99 


± 0.26 


^ 07 -1-0 11 
O.U / It U.ll 


Pt- An 


Q/l /I 8Q9 


8.80 


± 0.18 


cc 01 -I- 1 O 
O.Ui It U.iU 


i\/rr>i\/r 

iViiJlVi 


Q/I/1R /I /1 8 
044D.440 


8.64 


± 0.29 


/I Q1 -I- 1 O 

4.y± It U.iU 


Pt- AO 


3446.816 


8.63 


± 0.17 


4.97 ± 0.10 


MDM 


3450.832 


8.19 


± 0.16 


4.53 ± 0.09 


MDM 


3451.840 


8.38 


± 0.17 


4.36 ± 0.09 


MDM 


3452.793 


8.66 


± 0.17 


4.27 ± 0.09 


MDM 


3459.848 


8.99 


± 0.18 


4.56 ± 0.09 


MDM 


3460.809 


8.71 


± 0.17 


4.58 ± 0.09 


MDM 


3461.812 


9.00 


± 0.18 


4.56 ± 0.09 


MDM 


3462.848 


9.00 


± 0.18 


4.63 ± 0.09 


MDM 


3463.425 


10.23 ± 0.34 


4.78 ± 0.10 


CrAO^ 


3464.426 


9.58 ± 0.32 


4.76 ± 0.10 


CrAO^ 
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Table 1 — Continued 



JD 


Fx (5100 A) 


H/3 A4861 


Data Set 


(-2450000) 


(10~^^ erg A"^) 


(10 erg s ^ cm "^j 




3465.746 


9.19 ± 0.18 


4.80 ± 0.10 


MDM 


3465.890 


9.93 ± 0.15 




MAGNUM 


3467.844 


9.19 ± 0.18 


4.70 ± 0.09 


MDM 


3469.355 


9.33 ± 0.31 


4.87 ± 0.10 


CrAO 


3469.785 


9.41 ± 0.19 


4.80 ± 0.10 


MDM 


3470.364 


9.34 ± 0.31 


5.05 ± 0.11 


CrAO 


3470.852 


9.54 ± 0.19 


5.00 ± 0.10 


MDM 


3471.848 


9.29 ± 0.19 


4.90 ± 0.10 


MDM 


3478.823 


9.60 ± 0.07 




MAGNUM 


3485.007 


9.70 ± 0.07 




MAGNUM 


3508.910 


9.91 ± 0.05 




MAGNUM 


3527.847 


9.99 ± 0.03 




MAGNUM 


3539.818 


10.61 ± 0.12 




MAGNUM 


3556.769 


10.46 ± 0.06 




MAGNUM 


3572.766 


8.44 ± 0.12 




MAGNUM 


3579.766 


9.35 ± 0.25 




MAGNUM 



'^Continuum point omitted from final light curve. See §3. 
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Table 2. Light Curve Statistics 







Sampling 




Mean 






Time 




Interval(days) 


Mean 


Fractional 






Series 


N 


(T) 


^median 


Flux^ 


Error 


F 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


5100 A 


25 


1.7 


1.0 


4.8 ±0.7 


0.06 


0.14 


1.87±0.19 




27 


1.6 


1.0 


8.4 ±0.5 


0.02 


0.05 


1.25 ±0.04 



^Same flux units as Table 1 for 5100 A continuum and H/3, respectively. 
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Table 3. Reverberation Results 



Parameter 


Value 


(1) 


(2) 


Tcent 


3.73 ± 0.75 days 


''"pecik 


3.41°:^ days 


criine(mean) 


1790 ± 3 km s"^ 


FWHM (mean) 


5143 ± 16 km s"^ 


crii„e(rms) 


1561 ± 55 km s^^ 


FWHM (rms) 


4141 ± 416 km s'^ 


Mbh 


(9.8 ±2.1) X IO^Mq 
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Fig. 1. — Mean and RMS spectrum of NGC 4593 from MDM observations. The solid line 
shows the spectrum with the narrow-hne components of H/3, [Olll] A4959, and [Olll] A5007 
removed. The dotted line shows where these narrow-line components contributed to the 
spectrum before they were removed. The large increase in rms flux shortward of 4800 A is 
due to variations in the broad He ll A4686 emission line. 
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T 1 1 1 1 r= — I 1 1 1 1 1 1 1 1 1 1 r 




I I I I I I I I I I I I I I I I I I I I 

3400 3450 3500 3550 

JD (-2450000) 

Fig. 2. — Light curve showing complete data sets from all three sources. The top panel 
shows the 5100 A continuum flux in units of 10~^^ erg s~^ cm~^ A"^, while the bottom is 
the H/3 A4861 line flux in units of 10~^^ erg s~^ cm~^. The crosses show the two continuum 
points that were later omitted from the final light curve (see §3). 
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Fig. 3. — Light curve showing subset of data constrained to the time frame of the MDM 
observations that was used for time series analysis. Observations within 0.5 day of each 

other have been averaged together. The top panel shows the 5100 A continuum in units of 
Qj-g g-i cni-2 A~^, while the bottom is the H/3 A 4861 line flux in units of 10~^^ erg 



— 1 —2 

s cm . 
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T(days) 

Fig. 4. — Cross correlation function (CCF), discrete correlation function (DCF), and auto 
correlation function (ACF) from time series analysis of the continuum and H/3 light curves 
of NGC 4593. 
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Fig. 5. — Cross correlation functions (CCF) from time series analysis of the red and blue 
sides of the broad B./3 emission line (solid line) and of the wings and the core of the H/3 
emission hne (dashed line) of NGC 4593. 



